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1 Introduction 



Matrix models with the coupling to external matrices plays an important role in the contem- 
porary mathematical and theoretical physics. Historically, the first model of such type was 
the Brezin-Gross (BG) model ]I[ of the unitary matrix U linearly coupled to an external 
matrix field A. But the real breakthrough in this field was caused by Kontsevich's papers || 
where the generating functional for the 2D topological gravity was proved to be the inte- 
gral over the Hermitian matrices X with the potential X 3 , which are linearly coupled to 
an external matrix A. Simultaneously, the Witten hypothesis that this generating func- 
tional is a r- function of the KdV hierarchy was proved @, ||. The generalized Kontsevich 
model (GKM) — the model with an arbitrary polynomial potential V(X) and coupling with 
an external field — turned out to be a r-function of the Kadomtsev-Petviashvili hierarchy ||. 

Then, the interest to matrix models with logarithmic potentials appeared. The first 
such model with the external field coupling was proposed in || (the authors named it the 
Kontsevich-Penner (KP) model) and was pushed forward in [|7|, where its equivalence 
to the Hermitian one-matrix model with an arbitrary nonsingular potential was proved. 
Underlying geometrical structure is the discretized moduli space (d.m.s.) construction ||. 
Later on, the exact relation was proved that connects this model in the d.m.s. times with 
two copies of the Kontsevich integral taken at different time sets |10|| . 

Both the Kontsevich and the KP models, as well as the BG model, admit an explicit 
solutions in the 1/X-expansion [Q, ||, ||. Such solutions arise from the loop equation (or 
the Virasoro algebra constraints), which are at most quadratic in fields. One can formulate 
the problem to find all external field matrix models that manifest this property. Another 
model of this kind was the so-called NBI matrix model of IIB superstrings with the potential 
XA + X -1 + (2?7 + 1) logX appeared O, [12j in the context of the (M)atrix string theory. 
This model includes the BG model as a particular case (rj = 0) |y| and away of this point, 
it can be reduced [TJJ], after the time changing, to the Kontsevich model. (In particular, this 
enables one to produce the answer for the NBI model in the moment technique as soon as 
the answer for the Kontsevich model is known.) Note that the proof of equivalence of these 
two models relies upon the coincidence of the Virasoro algebras. 

The last model, which completes the list of matrix models with the loop equations 
quadratic in fields and which can be therefore solved in the 1/X-expansion framework is 
the two-logarithm (2-log) model with the potential XA + a log(l — X) + (3 log(l + X). This 



model turns out to be closely related to the exactly solvable Kazakov-Migdal models ITS 



and it was thoroughly investigated in the case of the unit matrix A, i.e., where it is reduced 
to the one-matrix model. Even in this case, this model manifests a rich phase structure Wfi 



In the present paper, we do not investigate all possible phases of the 2-log model and 
rather confine our consideration to the Kontsevich phase only, in which the expansion over 
traces of negative powers of the matrix A makes sense. First, we solve this model in the 
leading order of the 1/X-expansion; then, we find the constraint equations (the Virasoro 
algebra) and prove that in the proper normalization, these equations are exactly equivalent 
to the constraint equations of the KP model 0. Possible applications of the 2-log model are 
discussed. 



1 



2 Matrix model with two logarithms 



We start with the following matrix integral, which appear, for instance, in the logarithmic 
Kazakov-Migdal model [TBI, IT5|]: 



Z =JdXe ~ N tr [XA+a lo s( 1 " x )+' 3 1 °g( 1 +^)] . (2.1) 

This integral is of the most general form, since, rescaling and shifting the fields X and A, 
we may change the logarithmic branch points; however, we cannot change the constants a 
and P, which are actual charges in the model (|2.1|) . 

The matrix integral ( |2.1| ) belongs to a class of generalized Kontsevich models (GKM) ||T7f . 
Such models with negative powers of the matrix X have been previously discussed in the 
context of c = 1 bosonic string theory JT8 . 

For the models of this type, the large N solution is known explicitly only in some special 
cases. The models with cubic potential for X ||T^[ and the combination of the logarithmic 



and quadratic potentials || |2(J were solved by a method based on the Schwinger-Dyson 
equations, developed first for the unitary matrix models with external field |], |2l| . The same 
technique, being applied to the integral Q2.1|), also allows one to find its large N asymptotic 
expansions in the closed form for arbitrary a and f3. 

The Schwinger-Dyson equations for (|2.1|) follow from the identity 



l) (> _l^V dX ^e- Wtr [ XA+alos ( 1 - x ) +/31og ( 1+ ^ =0. (2.2) 



\N 3 dA jk dk u N JK 11 J J dXij 
Written in terms of the eigenvalues, these N equations read 




It is convenient to set 



Z{\)=0. 
(2.3) 

W(X^^lo g Z. (2.4) 



We also introduce the eigenvalue density of the matrix A: 

P(*) = JjY,8(x-* t ). (2-5) 

i 

The density obeys the normalization condition 

fdxp(x) = 1 (2.6) 



and in the large N limit becomes a smooth function. 

A simple power counting shows that the derivative of W(Xi) in the first term on the left 
hand side of equation fl2.3p is suppressed by the factor 1/N and can be omitted at N = oo. 



The remaining terms are rewritten as follows 

r2< \ [ j„ „/„ \ W(y)-W(x 



- xW 2 {x)-[dyp{y)y—^l VL±*1 + ( a + p _ 2)W(x) + {(3 - a) + x = 0, (2.7) 

J v — x 
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where Aj is replaced by x. Equation ( [2.7D can be simplified by the substitution 

W(x) = xW (x) - a + & - 1 . (2.8) 



After some transformations, using the normalization condition (|2.6|) , we obtain 

W\x) + x/dy p(y) ^ " = x 2 + (/3 - «)* + fct^ . (2.9) 

j y — x 4 

The nonlinear integral equation ( |2.9| ) can be solved with the help of the anzatz 

jy W = /W + ffc^)M^M, (2 . 10) 
2 •/ /(y) y-x 



where /(x) is an unknown function to be determined by substituting ( |2.10| ) into Eq. (|2] 
The asymptotic behaviors of W(x) and f(x) as x — > oo follow from Eq. ( |2.9|) : W(a;) ~ 
x + (/3 — a — l)/2, and the analytic solution with minimal set of singularities is merely 

f(x) = Vax 2 + bx + c. (2.11) 

Let us introduce the moments of the external field 

Hm* J °=Jwf dx - (2 - 12) 



The parameters a, 6, and c are unambiguously determined from Eq. ( |2.9| ). We find that 

c = ((3 + a - l) 2 /4, 
and a and 6 are implicitly defined by the following two constraints: 



i + \h = 4= , 

2 V a 

v^a«/o = /5 — a 

a 



(2.13) 



or, in terms of the eigenvalues, 



j 

So, we have 



i + ±y ! = ± 

a T7 ^ i ^ = /5 — a — - 

^ T Jatf + bXi + c a 



(2.14) 



Vax 2 + fcr + c 1/- /(y) - f(x) a + (3 - 1 

W{X) = X + 2J dyp{y) f(y)(y-x) + " (2 ' 15) 



Then, integrating fl2.15| ) w.r.t. x and checking that the stationary conditions w.r.t. the vari- 
ables a and b take place, we find the answer for the integral in the large N limit, 



logZ = N 2 ((3-a) 



1 1 

- log(6 2 — 4ac) — - log a 



+ N 2 {(3-a) 



^loga-^log(6 2 -4ac) 



+ y/c arctanh 



+ nj: 



2y/ca b 
~~b 2a 

a + (3-1 



+ N 2 



b 2 c 2c c. ,, 2 4 . 
■z-j ~ — + -7= + - log (b 2 - Aac) 
8cr 2a Ja 2 



\fc arctanh 



+ arctanh 



., logAi , 
2 y/a 

C+Xj b/(2y/c) \ 

f(\) ) 

'aX i X j + (X i + X j )b/2 + c 



^(j3-a) log (y/EXi + 7^ + f( X i 



E 



log(Ai - Xj) 



One can verify directly that 



<9 <9 

— \ogZ = — logZ = 



and ^- logZ = W(A{), as far as Eq. ( g3§ hold. 



(2.16) 



(2.17) 



3 A large N limit comparison 

Let us establish a relation between the constraint equations of the 2-log model and KP 
model ||. It is convenient to introduce new charges (parameters) instead of a and f3 



7 = (/?-a)/2, 



(a + (3-l)/2 = Vc, 



(3.1) 



and new variables 



b = b/a, c = c/a. 



(3.2) 



Shifting all eigenvalues Aj by the same constant £, we can rewrite the 2-log constraint equa- 
tions as follows: 



y/Z 2N 2 ? 



Af+(2£ + &) Xi+(e + & + 5 



A, + e 



(3.3) 



A 2 + ( 2f + b) Xi + (e + ib + c 



27-6, 



where "±" depends on the branch of the square root. Then, we make the following time 
change: 

tr^ = tr^ ± (- 2<p + 2N8 nA - mj) , (3.4) 
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where the role of "±" is the same. Making the presented time change and connecting the 
variables of the two models 



2£ + b = 46, 
£ 2 + fb + c = 4c, 
where 6 and c are already the KP variables, we have 

1 ^ 1 



(3.5) 



2b±—Y 

N ; 



rjf + Abrji + 4c 



0. 



(3.6) 



± — E 



+ c - 3b 2 = 7 - <f, 



i.e., exactly the constraint equations of the KP model with 7 — if = a. 0. Here £ is an 
arbitrary parameter. Using the original parameters a, (3, and a KP (a KP is the parameter a 
of the KP model, and a KP + 1/2 = a in the notation of 0), we see that a KP — (3 — 1. 

Naively, the parameter (3 is more preferred than a for some reason. Indeed, they play 
equal roles. The obvious symmetry of the 2-log matrix integral is encoded in the transfor- 
mations Aj — > — Aj (z = 1, N) and a /3. Under such a symmetry, 7 — > —7, ip^-if, and 

«kp = 7 - V 2 — ► «kp = -7 - V = a - 1- 

Let us recall the answer in the large N limit for the KP model 0. Substituting a = j — <f, 
we have 



log^K 



N 2 



7 - V 



I) log(6 2 -c)-^ 2 c-(7-^)c+^ + 3(7-^)6 2 + ^ 4 



'77. 



+ 6^ + c + (7 - pjlog 77i + 2b + + bm + c + 



4 E log I 



, ■ 2 fa + %0 + c + V "J + b ^ + c V "4" + + c 



(3.7) 



Here we compare the large N limit answer for the 2-log model with ( \5.7\ ). Further all 
equalities hold up to pure complex constant and irrelevant factors which can polynomially 
depend only on the parameters a and (3 (the polynomial of no more than second degree) of 
the 2-log model. Obviously such additional terms cannot influence the critical behavior of 
the model. Making the eigenvalue shift by £ and denoting d = b 2 — c, we obtain 

N 7 ' ' Ar2 ~ ' ^+^logd + 2^1og(6 + 2Vg)-6 



logZ 



log d + iV 2 7 



+ N 2 



2d + 2ipVc 



\)lo g d- if (ip + I) logc] + iVEflogfl + f + ^/l+J + f) 



' 46 4c 6 
Ai A 2 2Vd 



log Ai + Aj 



E/ 



A, 
1 1 

Aj Aj 



+ \/A? + 46Aj + 4c - A, 



(3.8) 
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where 

/ (x, y) = log (- + b (x + y) + 2cxy + -y/\ + Abx + Acx 2 yl + Aby + Acy 2 \ . (3.9) 
After some tedious algebra (similar to the one in f[4||), we obtain 

log Z = log Z KP + N ]T j (7 - if - I) log (^) + A; - y I + 2<^ 2 log p. (3.10) 

The difference between log Z and log Z KP depends only on some normalization factors in 
the large N limit. As is worth noting, these factors differ from the original normalization 
factors of the two models, which can be obtained by the early developed scheme 0]. We 
show that the appeared normalization factors are indeed natural. 

Let us investigate the Kontsevich regime of the two models (A — > 00 and rj — > 00). Then, 
for KP model we have (up to a constant) 

Z KP =jDXe Ntr [v X-4 +a lo g X] = e f tr„* (det v)a N 7^ tftr[-£+«lag(l+f )] ^ f ^ ^aiV 

There are two stationary points, Xo = ±1 + y/A, in the Kontsevich regime for the 2-log 
model (Y is the new variable). Choosing X = — 1 + Y/A for definiteness (another stationary 
point gives the same answer after symmetry A — > —A), we obtain 

Z = (det A )NW-l) e mrAj DYe -mr[Y + alo g (2-r) +(3logY] „ ( det A ) AT(/3-l) e iVtr A_ 

This is nothing but our normalizing factors. 



4 Constraint equations 



d 

Let us make the eigenvalue shift in the master equation of the 2-log model (di = 

-^ 2 (\ + Od^-^T,^l(9 ] ~d l ) + a+ ^~ 2 d l + (3-a + X l + ^Z = 0. (4.1) 
Using our normalizing factor 



(4.2) 



and pushing it through derivatives, we replace 

Then, we obtain master equation for the normalized partition function 

1 ft + 9l4EM(V*)-^ + ^s, 



N 2 



n 2 fr, a, - \ 



2^-1) ^-1)2 ^-1 



ATA,- 



A 



A; 



Z = 0. 



(4.3) 



(4.4) 
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Let us introduce the times of the 2- log model 

t n = -Y-!-- (4.5) 

Then, the constraint equations for Z({t n }) are obtained after some tedious algebra which 
we omit here. Collecting all coefficients to the term l/(AfiV 2 ), we obtain 

L k Z({t n }) = 0, k > -1, (4.6) 



where 



L k = V k+1 + £V k + £N(a + (5 - 1) ((1 - 4,o - - " ^m) 

+ C4,-iiV(/5-l)(ti-2iV), ' (4.7) 



and 



OO Q fc-1 Q Q Q 

V k = -T, mt ^ f E^T Jf f N(a - (3 4,o - 8k,-i) W 

m=l Ui ra J rk m=l m Ul k—m Ul k 

+ 27V(1 - 5^)-— + hS^-— + N 2 a((3 - l)4,o • (4.8) 

Ctfc+l <Jtk+l 

Here, the derivatives over to and t-i are fictitious and are used only for compactifying the 
presentation. 

For k,l > —1, L k satisfy the algebra 

[L k , Li] =(l- k)(L k+l+1 + £L k+l ). (4.9) 

Zero shift (£ = 0) results in the Virasoro algebra where the L_i generator is absent, 

[V k ,Vi] = (I - k)V k+b k,l>0. (4.10) 

We can also obtain the Virasoro algebra from the general algebra with nonzero shift by the 
replacement 



OO 



^ = EVr L ^' k>-l, (4.11) 

s=0 S 

which is singular at £ = 0. Performing the replacement and using the relations a KP = (3 — 1 
and if — — (a + /3 — l)/2, we obtain 

00 d k ~ 1 d d d d 

U = - £ mt m -f— - E j£— + 2Na° + 2N-^- 
m= i ot m+k m=1 Ol m Ot k - m Ot k Ot k+ i 

00 1 (9 d 
+ Vl ( il _ 2JV _^)( JVaKp + _|_). (4.12) 



After the time changing 



N N 
t n = t n - 2y? ( _^ ra + 2N5 n> i - y5 n , 2 , 



(4.13) 



where 



(4.14) 



are the times of the KP model, we obtain 

00 ~ d fe_1 d d d d 

C k = -J2 ™t m -~ E or- -sr— + S^kp^ + 

m=1 Otm+k m= i Ol m Olk-m Oik Olk+2 

- 2Na KP (5 kfi + 4,-i)^r ~ ^ 2 «KP*fc,o + 7T— + Na^hh,-! ■ (4.15) 

This is exactly the Virasoro algebra that appears in the KP model with the normalizing 
factor 

(4.16) 



Lb 



Indeed, we can perform the same operation for the KP model. First, we write the master 

d 

equation for the normalized partition function (di = — — ) 



dr]i 

1 °j ~ °i Vi a 2a KP Q i a KP ^ 1 a£ p 



——fp - \ ~* -1 — — f) — KP p> _i_ Kp \ ^ _ 



0. 



(4.17) 



Then, using the KP model times t n and collecting all coefficients to the term l/(rj^N 2 ), we 
obtain 

C k Z KP = 0, k > -1. (4.18) 

Therefore, we have proven the equivalence between the 2-log and KP models. 

Now, we write the explicit relation between the normalized partition functions of the two 
models 



'2- log 



-t4 \;a,p] = C(a,P)e^- 1)N2 e N2 W-^Z KP [f n (^), 
n \ n 



(),. 



(4.19) 



where 



i n (t,<p),a 



Z 2 -io S [A; a, f3] 



(4.20) 



n,{(^ Kp ef^ 2 }' 



a KP = (3 — 1 and C(a, (5) is some constant depending on the parameters a and f3. 

Note that we use here unshifted initial field A and explicitly show the dependence on the 
arbitrary parameter £ by the following reason. For the unshifted A-field, the Virasoro algebra 
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for the 2-log model does not possess the L_i generator. So, a question arises how we can 
obtain the £_i generator of the KP model when passing to the KP model. The reason is that 
after the time changing ( |4. , the KP times t n become ^-dependent. Differentiating Q4.19 ) 
over £ and using the relation 

^ = (n + 1)4+1 - N8 nA , (4.21) 
we obtain one more equation for Z KP , 

C-iZ KP = 0, 

where £_i is just the generator of the KP Virasoro algebra. 

5 Higher genus expressions 

Let us recall the genus expansion for the KP model ||, 

oo 

logZ KP = £iV 2 - 2 ^, (5.1) 

9=0 

where 

r-i / a a\ a \ M ai . . . M a J p 1 . . . J Q r r n\ 

f 9 = £ («i fa.--- fa\<xfa) g — — Tf3 Z — -> 9>h (5-2) 

and 

F^-^MMiJ^ 4 ), (7 = 1. (5.3) 
The moments were defined as follows (k > 0) 



1 w 1 

Mk = N § (?7i - x+) fc+1 /2 (r/i-x^) 1 ^ ~ ^ 



1 A 1 



Jfc ^ S fa " ^+) 1/2 (7ft - x_) fc +V2 ^ 



^^(A,-i/ + ) fc + 1 /2(A i - 2/ _)i/ 



2 



1 A 1 



(5.4) 



where x± are the endpoints of the cut for the one-cut solution and d = x + — x_ . In our 
notation, 

x± = -2b ± VAb 2 - 2c . (5.5) 
Let us introduce the moments for the 2-log model (k > 0) 

1 A 1 



(5.6) 
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where 



y± 



b -±J b - 

2 V 4 



c . 



(5.7) 



We are interested in the relation between the moments of the two models for k > (for k = 
the relation is given by constraint equations 



and (|3.6| )). After making the eigenvalue 
shift (y± =x± + £) and performing the time changing, we obtain (k > 1) 



(_l)fc+i 
M k = N k + 2 V \ [ 

y+ v- 

J k = K k + 2cp 1/2 fc+1/2 . 



y+ y- 



So, for the 2-log model, we have 



(5.? 



where 



and 



logZ = ^ N 2 ~ 29 Fp° e , 

9=0 



53 (oti . . . a s ; A • • • Pi\aPi) g n(N ai + 2<p 



aj>l,/3j>l 
I 



i=l 



-1 



ai+1/2 1/2 



X 



i=l 



n 



1/2 fli+1/2 

y+ y- 



3/2 1/2 

y+ y- 



/ 1 \P s_l 

x^i + 2y 1/2 3/2 j (y+-y-) 7 } , 



£ > 1, 



(5.9) 



(5.10) 



24 v y+ y- y v y+ y- y 

Therefore, expression (|2.16| ) for genus zero, taking into account the normalizing fac- 
tor ( f4.2| ), and expressions ( |5.10| ), ( |5.11| ), completely determine the partition function of the 
model (2.1) at all genera. 



6 Determinant formulas 



The exact determinant formulas in our model can be easily found using the Itzykson-Zuber- 
Mehta technique for the integration over angular variables in multi-matrix models. The 
partition function can be expressed as follows 

Z = (2n)^ /n [dxiil - Xi y aN (l + x t )-^e ~ NX ^} ^ , (6.1) 
where A (a;) = IIj> 7 ( 

the Van der Monde determinant and 6\ 2 are some integra- 
tion limits. We know that in the large N limit, the difference between partition functions 
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calculated in various integration limits is exponentially small and does not affect the 1/N 
perturbative expansion. So, we investigate several cases. 

(i). For 61 = — 1 and 9 2 = 1, we use the following integral representation (a, b > 0) 

l 

y c /x(l-x) a - 1 (l + x) fe - 1 e- cx = 2 a+b - 1 e- c B(a,b)<5>(a,a + b;2c), (6.2) 
-l 

where $(a, c; z) = i-F\(a, c; z) is the confluent hypergeometric function and B(a,b) is the 
beta-function. 

Then, in the domain a, [3 < 1/N, we obtain 

Z = {2 7l )^2-^ N2+N( - N ^/ 2 H^Bi-aN + 1, -f3N + *)} 



e -ATtr A 



x "TTvT , de L ||$(-aiV + 1, -(a + /3)iV + j + 1; 2iVA,)|| . (6.3) 
(ii). For 9 1 = 1 and # 2 = oo, we use the relation (a, c > 0) 

-s) 0-1 ^ + x) b - 1 e~ cx = (-l)°- 1 2° +6 - 1 e- c r(a)*(a,a + 6;2c), (6.4) 



i 

where 

tt(a, c; z) = -IS— j- $(a, c; *) + ^f^- $(a - c + 1, 2 - c; z) (6.5) 
r(a — c+1) r(a) 

is the other confluent hypergeometric function and T(a) is the gamma-function. 
Then, in the domain where a < 1/N, (3 is unrestricted, and A, > 0, we have 

Z = ( 2vr )^(-i)-^ 2 2 -(-+^ 2 +^+ 1 )/ 2 T N (-aN + l) 

—Ntv X 

x ~7\7rT ^ de k, I ^(-^ + !> "(« + Z 3 )^ + i + 1; 2A ^)I I • ( 6 - 6 ) 
A(A) i<*,:?<2V 



This answer covers more general domain of the parameters a and /3 than 

(iii). If a = we get the simplest answer setting Q\ = — 1 and 9 2 = 00. In the domain 
f3 < 1/N and Aj > 0, we obtain 



a=0 



(2vr)^ J] {r(-/?iV + ^} (det \)(P-^ e NtrA ^ (6J) 

i 

which is the unshifted normalizing factor up to a constant. 
7 String susceptibilities 



Let us calculate the string susceptibility w.r.t. 7 and (p for (|2.16|) . By virtue of Eq. (|2.17|) , 

d d 

— log Z = — — log Z and the same holds true for (p. Furthermore, an amazing fact is that 

cry cry 
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the expressions obtained are themselves stationary w.r.t. differentiation over a and b. This 
means that the total second derivatives in 7 and ip coincide with the corresponding partial 
derivatives, so we have 

1 d 2 

X 1 = J^2~^2 lo S Z = lo s( &2 -4ac) -2 log a, 

Id 2 2Jac b + 2Jac 

* 2 = ~at2 1~T~ § = -2arctanh— — = - log- — — = , (7.1) 
N z d^dip b b — 2a/oc 

1 d 2 

X3 = j^-^ lo s z = lo s( fo2 ~ 4ac ) + 6 - 

Recalling the string susceptibility of the KP model in the KP variables b and c || 

Xkp = log(6 2 - c) (7.2) 
and using the relations ( |3.5| ), we obtain 

Xkp = Xi- (7-3) 



8 Conclusion 



This paper concludes the series of papers || |10|, [12], devoted to styduing the external field 
matrix problems with logarithmic potentials. We see that, at least in the 1/iV-expansion in 
terms of the corresponding moments, all these models can be reduced either to the Kontsevich 
model or to the Hermitian one-matrix model with an arbitrary potential. Here, the question 
arises whether this can be derived directly within the r-function framework || . The related 
question is which reductions of the Kadomtsev-Petviashvili hierarchy correspond to the NBI 
and 2-log models. 

One can always say the origin of the logarithmic terms is due to additional degrees of 
freedom that were integrated out. Matrix integral (|2.1|) can be represented as the 0(a, 0)- 



type |24j matrix integral 



-Jdxu 

J i=l 



dVidWi _Q d^jd^je L 31 \ , (8.1) 



where the sum over repeated indices is implied and we assume the matrix fields $ and ^ are 
Grassmann even. Action ( |8.1| ) is of a nonlinear sigma model type with free matrix fields $ 
and ^ dwelling on the manifold A — + = 0. 
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